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ABSTRACT 

We measure the autocorrelation function, £, of galaxies in the IRAS Point Source 
Catalogue galaxy redshift (PSCz) survey and investigate its dependence on the far- 
infrared colour and absolute luminosity of the galaxies. We find that the PSCz survey 
correlation function can be modelled out to a scale of 10 /i _1 Mpc as a power law of 
slope 1.30 ± 0.04 and correlation length 4.77 ± 0.20. At a scale of 75 ft^Mpc we find 
the value of J 3 to be 1500 ± 400. 

We also find that galaxies with higher lOO/im/60/im flux ratio, corresponding to 
cooler dust temperatures, are more strongly clustered than warmer galaxies. Splitting 
the survey into three colour subsamples, we find that, between 1 and 10 /i _1 Mpc, the 
ratio of £ is factor of 1.5 higher for the cooler galaxies compared to the hotter galaxies. 
This is consistent with the suggestion that hotter galaxies have higher star-formation 
rates, and correspond to later-type galaxies which are less clustered than earlier types. 

Using volume limited sub-samples, we find a weak variation of £ as a function 
of absolute luminosity, in the sense that more luminous galaxies are less clustered 
than fainter galaxies. The trend is consistent with the colour dependence of £ and the 
observed colour-luminosity correlation, but the large uncertainties mean that it has a 
low statistical significance. 

Key words: galaxies: distances and redshifts - infrared: galaxies - large-scale struc- 
ture of Universe 



1 INTRODUCTION 

The relation between the distribution of galaxies and the 
distribution of mass is now one of the most important prob- 
lems in large-scale structure. Empirically, it has been found 
that different galaxy types have different clustering ampli- 
tudes and hence that galaxies are biased tracers of the mass 
distribution. Evidently, we need to understand the physical 
mechanisms responsible for these biases if we are to establish 
the connection between galaxy tracers and the mass. This 
problem can be understood in terms of two related ques- 
tions: (i) How does the clustering of a galaxy sample de- 
pend on the properties of the galaxies? and (ii) How do the 
properties of galaxies depend on their local environment? 

Many ideas have been put forward as to how the en- 
vironment may modify galaxy properties, from schematic 
ideas concerning feedback in the formation process (White 
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& Rees 1978; Dekel & Rees 1987) to hydrodynamical simula- 
tions of galaxy formation incorporating cooling and dissipa- 
tion (e.g. Katz & Gunn 1991; Navarro & White 1993; Pearce 
et al. 1999). Over the next few years we can expect ma- 
jor advances in hydrodynamic simulations employing paral- 
lel computers and there is a real prospect of understanding 
the relationship between galaxy morphologies and environ- 
ment within the context of specific theories of cosmological 
fluctuations. We can hope that such numerical simulations 
will eventually be able to provide detailed predictions of the 
spatial distribution of galaxies as a function of their stellar 
content, rotation velocity, star- formation rate, and morpho- 
logical type. 

The simplest bias schemes are those based on the statis- 
tics of high density peaks (Kaiser 1984, Mo & White 1996). 
Since the optical magnitude of a galaxy is correlated to the 
depth of the potential via the Tully Fisher relation for spirals 
and D n — a relation for ellipticals, this would suggest that 
the brightest galaxies will be more strongly correlated than 
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fainter galaxies. Realistic galaxy formation models rely on 
feedback mechanisms to limit star-formation, either acting 
internally within each galaxy or involving interactions with 
other galaxies and the intergalactic medium. So galaxies in 
regions of high local galaxy density are likely to have a re- 
duced star- formation rate (SFR). Indeed this is confirmed 
by the observation that later-type optical galaxies, which 
have higher SFRs, have a much lower clustering amplitude 
than earlier galaxy types (e.g. Rosenberg, Salzer & Moody 
1994; Loveday et al. 1995; Loveday, Tresse & Maddox 1999). 

Galaxies with a high SFR tend to be bright in the far 
infra-red (FIR) because of the thermal emission from dust 
heated by young stars. Thus the observations that IRAS 
galaxies (selected on their 60/xm flux) tend to avoid rich 
galaxy clusters, and have a lower clustering amplitude are 
consistent with this general picture. In this paper we divide 
the PSCz sample of IRAS galaxies into subsamples based on 
their FIR colour. Since the galaxies with warmer FIR colours 
are those with higher SFR, we can further test the de- 
pendence on star-formation rate. Mann, Saunders & Taylor 
(1996) selected warm and cool sub-samples of IRA S galaxies 
from the QDOT survey based on their 60/im/100/im flux ra- 
tio, and found marginal evidence that warmer galaxies are 
more strongly clustered than cooler galaxies, the opposite 
to what is expected from the morphology density relation. 
However they did not consider the result very significant 
compared to the expected cosmic variance. 

Using optical galaxy samples, a number of authors have 
found indications of changes of the clustering amplitude with 
galaxy luminosity, (e.g. Valls-Gabaud, Alimi & Blanchard 
1989; Park et al. 1994; Moore et al. 1994; Loveday et al. 
1995; Benoist et al. 1996) but the samples are small and 
the observed amplitude shifts do not have a high statis- 
tical significance. The expected luminosity dependence for 
IRAS galaxy samples is less clear, since the FIR luminosity 
of a galaxy is likely to depend on both the mass of its dark 
halo and its specific SFR. If the FIR luminosity is corre- 
lated to the halo mass, then brighter galaxies should have a 
higher clustering amplitude. On the other hand, a high SFR 
will make a galaxy more luminous, and galaxies with a high 
star-formation rate tend to have weaker clustering. It is not 
easy to predict which of these is the dominant effect. Ob- 
servationally Szapudi et al. (2000) and Beisbart & Kerscher 
(2000) have used mark correlation functions with volume 
limited sub-samples to examine the luminosity dependence 
of clustering in the PSCz survey. Over the narrow range of 
luminosities they consider they concluded that there is no 
significant luminosity dependence. 

In this paper we investigate clustering amplitude for dif- 
ferent sub-samples of galaxies from the PSCz survey. We be- 
gin by summarizing the relevant aspects of the PSCz survey 
and describe various sub-samples that we will study here. In 
Section || we show our methods and in Section ^| we give our 
results for the whole sample and consider how the amplitude 
of the redshift space autocorrelation function depends on the 
colour and luminosity selection of the galaxies. Finally, we 
discuss our conclusions in Section H. 



2 THE PSC REDSHIFT SURVEY 



2.1 Description 

Galaxies for the PSCz were selected from the IRAS Point 
Source Catalogue (Beichman et al. 1988) by positional, iden- 
tificational flux and colour criteria designed to accept almost 
all galaxies while keeping contamination by Galactic sources 
to an acceptable level. The significant differences compared 
with the QDOT survey (Lawrence et al. 1999), were that the 
colour criteria were relaxed to avoid excluding previously 
uncatalogued galaxies with unusual (especially cool) colours 
and the sky coverage is increased to 84 per cent by including 
areas of high cirrus contamination. These changes increase 
the contamination of the galaxy catalogue from local galac- 
tic sources (stars, planetary nebulae, cirrus etc) but ensure 
the sample has a higher completeness. The contaminating 
Galactic sources were excluded by a combination of IRAS 
and optical properties (from sky survey plates), and where 
necessary spectroscopy. Redshifts were taken from all avail- 
able published or unpublished sources; principally Huchra's 
ZCAT, the LEDA database, the 1.2Jy and QDOT surveys. 
Also 4500 new redshifts were taken on the Isaac Newton 
Telescope, the Anglo- Australian Telescope, the Cerro-Tololo 
1.5-m and Cananea 2.1-m telescopes. In total the catalogue 
contains 14677 redshifts and the median redshift is ~ 0.028. 
The uniformity and completeness of the catalogue are dis- 
cussed by Saunders et al. (2000). 

The redshift distribution of the galaxies is shown in Fig- 
ure |l|. It can be seen that the peak of the distribution is near 
2 = 0.02, and there is a tail which extends out to beyond 
z — 0.1. This tail means that there are over 1200 galaxies in 
the survey with a redshift > 0.08. 

For our present analysis, only galaxies with a redshift > 
0.004, corresponding to a recession velocity of 1200 km s _1 
were selected, so that local effects and the velocity uncer- 
tainties of around 120 km s _1 were less significant. 



2.2 Sub-samples 

2.2.1 Colour sub-samples 

In order to investigate the dependence of clustering on the 
temperature of the galaxies, we used the addscan 60/im and 
100/im fluxes (/6o and /ioo) to sub-divide the parent red- 
shift catalogue into three sub-samples: hot galaxies, with 
/loo/ieo < 1.7 ; warm galaxies, with 1.7 < /loo/jeo < 2.3; 
and cold galaxies, with fioo/fm > 2.3. These boundaries 
correspond to blackbody temperatures of around 31 and 
28K, and were chosen to give roughly equal numbers in 
each sub-sample. The actual numbers of galaxies in the 
sub-samples were 4452, 4388 and 4107 respectively. The 
mean /100//6O colour ratios for the hot, warm and cold sub- 
samples were 1.31, 1.98 and 2.91 respectively, corresponding 
to black-body temperatures of about 34K, 30K and 26. 5K 
respectively. The mean for the whole catalogue is 2.05, cor- 
responding to a black-body temperature of around 29. 5K. 

The redshift distributions for the sub-samples are plot- 
ted in Figure ^](a) and it can be seen that the cooler samples 
tend to peak at a slightly lower redshift than the hotter sam- 
ples. This reflects the correlation between colour and abso- 
lute luminosity as plotted in Figure ^: cooler galaxies tend to 
be fainter, and so are only seen nearby. Nevertheless, these 
differences in redshift distributions are smaller than the dif- 
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Figure 1. The redshift distributions for galaxy samples, (a) 
colour selected sub-samples listed in Table (b) absolute lu- 
minosity selected sub-samples listed in Table B. 

ference for luminosity selected sub-samples, as discussed in 
the next section. 

The angular distribution of these sub-samples are plot- 
ted in Figure ^. It can be seen that these samples show no 
obvious gradients as a function of position on the sky, and 
also that the cooler sample appears to be more strongly clus- 
tered. Figure |^ shows the projection of the galaxies on to a 
plane along the celestial equator. Again the cooler sample 
appears more clustered than the warmer samples. Though 
apparently quite significant, these visual impressions should 
be treated with caution. The cooler sample is shallower than 
the warmer samples, so the angular clustering would be 
stronger, even if the samples had the same spatial cluster- 
ing. Also the cool sample has a higher space density than 
the warmer samples, and this may enhance the visual ap- 
pearance of clustering. 

2.2.2 Absolute luminosity sub-samples 

The parent sample was also divided into sub-samples of ab- 
solute luminosity at 60/im: faint galaxies, with log 10 (Z/6o) < 





1 1 1 1 1 




III 












- ' ■ 






_f o 


•4SS3 

- 






o 




lit 












- 










00 









-0.5 0.5 1 



lo g( f ioo/ f 6o) 

Figure 2. Log plot of absolute luminosity at 60fim against the 
100fj,m/60fim colour. The vertical lines show the boundaries used 
to select our colour sub-samples, and the horizontal line shows the 
boundary between the luminosity samples. 

9.6; and bright galaxies with log 10 (Lgo) > 9.6, where Lqo = 
4Tvd 2 feo and d is the luminosity distance to the galaxy (we 
assumed Q — 1 and A = to determine the geometry, but at 
the distances considered here, cosmolqgical effects are neg- 
ligible.) The units of Leo are LQh~ 2 Nl Note that the usual 
definition of an ultra-luminous IRAS galaxy corresponds to 
log 10 (L6o) ^ 11.4 (Soifer et al 1987), and so the lower limit 
of our brighter sample is a factor of ~ 60 fainter than ultra- 
luminous galaxy samples. 

The bright sample has a higher mean redshift, and 
hence covers a much larger volume, but the galaxies are 
a more dilute sampling of the density field and so a larger 
number is needed to give a similar number of small separa- 
tion pairs. Hence the faint and bright samples were chosen 
to have 4719 and 8228 galaxies respectively. 

The redshift distribution for these sub-samples is shown 
in Figure [l](b). Since the parent sample is flux limited, the 
low luminosity galaxies have a sharply defined upper red- 
shift limit corresponding to the distance where the observed 
flux equals the flux limit. For the more luminous galaxies, 
the clustering signal is dominated by the volume at higher 
redshifts. This means that the clustering measurements for 
the two samples come from almost independent volumes, 
and so cosmic variance introduces a large uncertainty in the 
comparison. The clustering measurements for colour selected 
sub-samples are less affected by cosmic variance because 
they cover very similar volumes, and so their comparison 
should be more reliable. 

2.2.3 Volume limited sub-samples 

The cosmic variance between the luminosity sub-samples 
can be partly overcome by using volume limited sub- 
samples. To do this, the parent catalogue was split into 
sub-samples limited by distance and also by absolute lu- 
minosity. Galaxies within a sphere of radius x, centred on 
us, are included only if their absolute luminosity (Leo) is 

t h is the Hubble constant in units of 100 km s — 1 Mpc - 1 
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Figure 3. Plots of the position of galaxies on the sky for the three colour sub-samples listed in Table |l| (a) hot, (b) warm, and (c) cold. 
The line shows the Galactic equator and the obscuration around this is masked out in the analysis. 
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Figure 4. The galaxy positions projected onto the celestial equator for the full survey and the three colour sub-samples listed in Table | 
(a) all, (b) hot, (c) warm, and (d) cold. 



greater than 4ivx 2 fmin, where fmin is the flux- limit of the 
survey or sample (0.6 Jy for the full PSCz survey). Without 
galaxy evolution these samples will be homogeneous. Unfor- 
tunately the numbers of galaxies in these samples, shown in 
Table are rather small. 

IRAS galaxies are known to show very rapid evolution 
(Saunders et al 1990), and this could introduce a radial gra- 
dient in space density in the volume-limited samples. Fig- 
ure |E| shows the galaxy number density as a function of ra- 
dius, normalised so that the maximum radius of the vol- 
ume is 1, for the samples with radius 100, 150, 200 and 
300 ft -1 Mpc. No significant gradients are apparent, and so 
we assume a constant density in our clustering estimates. 
These sub-samples allow another test for a luminosity trend 
in the clustering as the larger volumes are limited to brighter 



galaxies than the smaller volumes. Note that the different 
volumes have a significant fraction of galaxies in common, 
so that the results for different volume scales are not com- 
pletely independent. 

2.3 PSCz mock catalogues 

Large N-body simulations from Cole et al. (1998) were used 
to create mock PSCz surveys. The simulations used the 
AP 3 M code of Couchman (1991) loaded with 192 3 parti- 
cles of mass 1.64 x 10 12 f2o/i~ 1 Mq in a box with a co-moving 
size of 345.6/i _1 Mpc. More details can be found in Cole et al. 
(1998). For this analysis three different cosmological mod- 
els have been considered. They are a spatially flat low den- 
sity universe model (ACDM) with Qq = 0.3, Qa = 0.7 and 
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r = 0.25, a spatially flat universe model (SC ' DMYq.2z) with 
tto = 1.0 and F = 0.25, and finally a spatially flat universe 
model (SCDM co be) with Qo = 1.0 and F = 0.5, normalised 
to match the amplitude of the COBE data. For pure CDM- 
type models the shape parameter, T — Qh. 

The galaxy catalogues were extracted from the numer- 
ical simulations by first identifying a population of objects 
with similar properties to the Local Group (LG). Then 
an observer is placed by considering the observational con- 
straints of peculiar velocity, local overdensity and shear. All 
points within a sphere of 120/i -1 Mpc radius around the ob- 
server are included, and the frame is rotated so that the 
motion of the observer is the same as the LG peculiar veloc- 
ity. Note that this sphere radius is less than the size of the 
PSCz survey. A luminosity was randomly assigned to each 
galaxy in the mock realizations to mimic the corresponding 
observational properties. The PSCz density is reproduced by 
using the PSCz selection function to reject points and also 
to randomly assign a flux to each point. Unsurveyed regions 
are rejected using the PSCz mask to give the same sky cover- 
age as the real catalogue. Each mass point is associated with 
a galaxy so the linear bias parameter, b — 1.0. Lastly the 
bias against early-type galaxies in the original IRAS survey 
is introduced by rejecting an appropriate fraction of galaxies 
in clusters the cores of clusters, which is meant to reproduce 
Dressler's morphology-density relation (Dressier 1980). The 
resulting mock surveys have approximately the same corre- 
lation function, redshift distribution and flux distribution as 
the real PSCz sample. 

We calculated the correlation function for the 10 re- 
alizations of each model, and used the standard deviation 
about the mean of so and 7 to estimate of the uncertain- 
ties in these measurements for the real sample. These cata- 
logues can also be split into the same bright and faint sub- 
samples as the real catalogue. Again, the standard devia- 
tion between the realizations gives an estimate of the un- 
certainty on the measurements for the real data, and so al- 
lows more reliable estimates of the significance level of any 
variations seen in the real data. Colours can be included 
in the mocks by considering Figure which we modelled 
as a linear relation between log 10 (/ioo//6o) and log 10 (Z/6o) 
with a Gaussian dispersion about the mean. Each galaxy 
can then be assigned a colour according to its luminosity, 
IogroCfW/eo) = Alog 10 (L 6 o) + B + C, where A = -0.117, 
B — 1.427 and C is a random number selected from a 
Gaussian distribution with variance of 0.2. Then colour sub- 
samples can be selected and tested in the same way as the 
real data. Since the luminosities were assigned at random 
to the mock galaxies, and the colours are based on the lu- 
minosities, the clustering in the mocks is independent of 
luminosity and colour. 



3 METHODS FOR ESTIMATING THE 

REDSHIFT SPACE AUTOCORRELATION 
FUNCTION £ 

There has been considerable interest in analysing the best 
way to estimate the autocorrelation function from redshift 
surveys (e.g. Landy & Szalay 1993; Vogeley & Szalay 1996; 
Hamilton 1997a,b). Reliable and unbiased estimates can be 
obtained by cross-correlating the real data with an artificial 



unclustered catalogue which has the same selection func- 
tions in angle on the sky and in redshift. We create this 
catalogue by generating random positions uniformly over 
the sky, and rejecting positions behind the mask which de- 
fines the boundaries of the real data. For each position we 
then assign a redshift from a distribution which matches the 
observed data. 

We investigated several methods to generate the red- 
shift distribution for the random catalogues: random shuf- 
fling of the positions on the sky relative to the redshifts; 
random shuffling followed by adding a further random ve- 
locity from a Gaussian distribution with a = 500, 1000 or 
1500 km s _1 ; and generating a random redshift distribution 
to match an analytic fit to the selection function (Saunders 
et al. 2000). For the full catalogue we found that random 
shuffling with 1000 km s _1 smoothing and the analytic fit 
gave indistinguishable results for { . Rather than calculating 
a best fit selection function for each sub-sample we simply 
generated appropriate sub-samples using the 1000 km s _1 
smoothed randomized distribution. For each random cata- 
logue we generated 12 times the number of random points 
as real galaxies for the whole sample, and 20 times as many 
for the sub-samples. 

3.1 Flux limited sample methods 

The estimate of £(s) we used (Landy & Szalay 1993) is given 
by 

(DP) - 2(DR) + (RR) 
C(S) = (RR) ' (1) 

where (DD) is the weighted sum over galaxy pairs with sep- 
aration s in redshift space; (RR) is the weighted sum of 
pairs with the same separation in the random catalogue; 
and (DR) is the weighted sum of galaxy-random pairs. In 
each of these sums, the pairs at separation s are weighted 
by a factor 

WlJ = (1 + 47m( Zi ) J 3 (a))(l + 4nn(z 3 )Ms)) ' (2) 

where n(Zi) is the space density of observed galaxies at the 
redshift z; for galaxy i, and Js(s) = J Q £(r)r dr (Efstathiou 
1988). 

To calculate £(s) we first calculated £(s) with a standard 
pair weighting of the sample and fitted a power law to the 
results for 1 < s < 10 h~ Mpc. We then recalculated the 
correlation function using J3 from this best fit power law, 
truncated so that the maximum J3 is 1500. It was found that 
only one iteration was needed to produce a stable result with 
a consistent J3 and £. The calculation of £(s) is relatively 
insensitive to the precise form of the J3 weighting employed. 
We then fit the final correlation function with a power law 
of the form: 

^=(£T- (3) 

We calculate the correlation function for mock cata- 
logues in the same way as for the real sample. 

3.2 Volume limited sub-sample methods 

For the volume-limited sub-samples a different technique 
was used to estimate the correlation functions. Following 
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the method of Croft et al. (1997), we used a maximum like- 
lihood approach that assumes that each distinct galaxy pair 
is an independent object. This technique uses the fact that 
1 + £(s) can be considered as the probability distribution 
for galaxy pairs as a function of separation. Thus a power 
law £ can be fitted directly to the galaxy pair distribution 
without the need to specify arbitrary bins of pair separa- 
tion. The resulting power law parameters tend to be more 
robust than fitting to a binned version of £, especially for 
small samples. The homogeneity of volume limited samples 
means that simple unit pair weighting yields the best results. 

If galaxy pairs were independent, the calculated likeli- 
hood contours would directly provide 1-a uncertainties in 
the best-fit parameters, equivalent to the Poisson error in 
the pair counts. Since galaxy pairs are correlated, the effec- 
tive number of pairs is reduced, and so the uncertainties are 
increased by a factor (1 + ^imJz). These corrected Poisson 
estimates still do not allow for any uncertainty introduced 
by the variance in mean density on the scale of sample size. 
We used a simple approach to include this effect, by con- 
sidering the variance in the mean density within the sample 
volume, var(<5n/n) = 3J3/R 3 , where 7? is the radius of the 
volume. This introduces an uncertainty in £, which we add 
in quadrature to the corrected Poisson estimates. 

For volumes with radius R < 100 h~ 1 Mpc we have 
also analysed the ten ACDM mock catalogues using the 
same method. The ACDM model was chosen because its 
clustering amplitude is the closest match to the observed 
amplitude. The standard deviation between results pro- 
vides an independent estimate of the uncertainties in sq 
for each sample. For the volumes we consider (with radius 
R > 50 fo _1 Mpc) the value of J3 is roughly constant, and 
we find that the two methods agree to ~10% if we choose 
J 3 = 1500. For volumes with radius R > 100 /i _1 Mpc the 
mock catalogues are too small, and so we have to rely on 
the analytical approximation to estimate the uncertainties. 



4 RESULTS 
4.1 Full sample 

Our estimate of £(s) for the full sample is shown in Figure ^| 
The power law fit gives s = 4.77 ± 0.20 and 7 = 1.30 ± 
0.04. The errors are estimated from the standard deviation 
between the ten realizations of the ACDM mock catalogues. 
We calculated the value of J3 and found it to be 1500 ± 400 
at a scale of 75 h~ 1 Mpc. The error on the J3 result is from 
the standard deviation of the mock realization results. 



4.2 Colour sub-samples 

Our estimates of the redshift-space correlation functions for 
the three colour sub-samples are shown in Figure (?]. The er- 
ror bar for each point shows the standard deviation between 
10 realizations of the ACDM mock catalogues. It is clear 
that over the range 1 < s < 10 h~ 1 Mpc the amplitude of 
£(s) is higher by about a factor of 1.5 for the cooler sub- 
sample. Fitting a power law to the data over the range of 
1 < s < 10 /i _1 Mpc gives parameters as shown in Table |l[ 
clearly confirming the increase in clustering amplitude for 
cooler galaxies. Figure ^ shows x 2 contours for the hot and 
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Figure 5. Galaxy number density plotted against normalised 
radius for the volume limited samples with radii labelled. 



Table 1. Power law fit parameters to the colour sub-samples 



Free parameter best fits: 
colour range 
/ioo//eo < 1.7 
1.7 < /100//6O < 2.3 
/100//6O > 2.3 

Fixed 7 at 1.30: 
colour range 
/100//6O < 1-7 
1.7 < /100//6O < 2.3 
/100//6O > 2.3 



so 

4.41 ± 0.29 
4.76 ± 0.30 
5.49 ± 0.29 



so 

4.43 ± 0.30 
4.64 ± 0.33 
5.38 ± 0.30 



7 

1.32 ± 0.06 
1.21 ± 0.06 
1.25 ± 0.05 



cold sub-sample, and suggests that trend is significant at 
about the 2a level. 

Since the best-fit so and 7 are correlated, we have also 
fitted the data using a fixed value for 7, set to the best fit 
value for the full sample, 7 = 1.30. The resulting so values 
for the three different colour sub-samples are listed in Ta- 
ble [l] and plotted as a function of mean colour in Figure ^. 
The trend of increasing so for warmer galaxies is clearly ap- 
parent at about the 3cr level. A simple linear regression to 
the points yields 

log 10 (s ) = 0.241og 10 (4^) + 0.61. (4) 
.ho 

Analyzing colour selected sub-samples from the mock 
catalogues shows no significant difference in amplitude be- 
tween the hot and cold sub-samples. This is as expected 
since we assigned the colours at random, but the measure- 
ments provide an important check that there are no subtle 
biases in our algorithms. 



4.3 Luminosity sub-samples 

The estimates for the two luminosity sub-samples are shown 
in Figure Juj. The error bars are derived from the scatter 
between ten realizations of the ACDM mock galaxy cat- 
alogues. It is clear that there is very little difference be- 
tween the two sub-samples. The power law fit parameters 
are shown in Table ^| These parameters show no significant 
difference when the errors are taken into account. 
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s /h Mpc 

Figure 6. Logarithmic plot of £ against s for the full PSCz sample. The dotted line is the best-fit power-law, fitted for 1 < s < 10 /i -1 Mpc. 
The inset shows a linear plot of § against s. The error bars are derived from the scatter between ten realizations of the AC DM mock 
galaxy catalogues. 



Table 2. Power law fit parameters to the luminosity sub-samples 



Free parameter best fits: 
luminosity range 
log{L 60 ) < 9.6 
log(Leo) > 9.6 

Fixed 7 at 1.30: 
luminosity range 
log{L 60 ) < 9-6 
log(L 6 o) > 9.6 



so 

4.96 ± 0.44 
4.80 ± 0.31 



so 

4.82 ± 0.54 
4.89 ± 0.35 



7 

1.21 ± 0.04 
1.36 ± 0.06 



For all mock catalogues we found no difference between 
the results for the faint and bright sub-samples within the 
errors. This is as expected since we assigned the luminosi- 
ties at random, but the measurements provide an important 
check that there are no subtle biases in our algorithms. 



Table 3. Maximum likelihood fits to the volume limited sub- 
samples with 7 fix ed at 1.30, and uncertainties estimated as de- 
scribed in Section |3.2j . 



Volume radius 


no. of 




so 


(fe^Mpc) 


galaxies 


(h- 1 Mpc) 


50 


1699 


9.41 


5.11 ±1.19 


75 


2151 


9.71 


4.91 ±0.65 


100 


2188 


9.94 


4.23 ±0.44 


125 


1961 


10.12 


4.90 ±0.38 


150 


1633 


10.28 


5.25 ±0.42 


175 


1399 


10.41 


3.86 ±0.53 


200 


1202 


10.53 


3.56 ±0.74 


225 


958 


10.64 


2.83 ±1.08 


250 


779 


10.74 


4.35 ±1.43 


275 


653 


10.83 


4.49 ±2.00 


300 


587 


10.91 


3.65 ±2.50 



4.4 Volume limited sub-samples 

Our results from the volume limited sub-samples are shown 
in Table ^ and displayed in Figure [ll]. Although there is 
clearly a hint of a trend of reducing amplitude with scale, it 
can be seen that the error bars are too large to say that this 
is a statistically significant relation. 

The minimum absolute luminosity in the volume lim- 
ited samples increases with the radius of the volume, so the 
mean luminosity also increases. Figure [t] explicitly shows 
the clustering amplitude as a function of the mean luminos- 
ity in the volume limited samples. 



5 CONCLUSIONS AND DISCUSSION 

We detect a colour dependence in the clustering of PSCz 
galaxies: the cooler galaxies are more strongly clustered than 
the hotter galaxies. Although the significance level of ampli- 
tude variations is not high, this does not mean that the 
effects are small: the clustering amplitude differs by a fac- 
tor of 1.5 between the warm and cool sub-samples. We do 
not detect any significant variation in clustering for luminos- 
ity selected subsamples. Volume limited sub-samples show a 
slight trend corresponding to a weaker clustering for brighter 
galaxies, but the uncertainties mean it is not significant. 
Our results are different to those of Mann, Saunders 
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Figure 7. Logarithmic plots of § against s for the colour sub-samples defined in Table |l[ The inset shows a linear plot of £ against s. 
The error bars are derived from the scatter between ten realizations of the ACDM mock galaxy catalogues. For clarity, error bars are 
plotted only for the warm sub-sample. 




s„ Figure 9. Clustering length so versus mean colour for 7 = 1.30. 

The filled points show the colour sub-samples as defined in Ta- 
Figure 8. Contours of \ 2 as a function of sq and 7 for thepower- ble |l], and the open point shows the full catalogue. The best fit line 

law fits to £ of the hot and cold sub-samples as in Table 0. is also shown and has the form log 10 (so) = 0.241og 10 ('y m )+0.61. 



& Taylor (1996) who found marginal evidence that warmer 
galaxies are more strongly clustered than cooler galaxies. 
The main reason for the difference is simply the factor of 6 
increase in number of galaxies in our sample, which dramat- 
ically reduces the uncertainties in our clustering estimates. 
The uncertainties on the earlier measurements encompass 
our results, so there is actually no discrepancy. 



The variation of clustering with FIR colour that we 
find in the PSCz is consistent with the variation with star- 
formation rate indicated by Ha and O [II] emission lines in 
the SAPM survey (Loveday et al. 1999). Loveday et al. split 
the SAPM survey into three samples with low, medium and 
high Hct equivalent widths (EW), and estimate so for the 
samples to be 8.7, 5.5 and 4.6 /i -1 Mpc respectively. A simi- 
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Figure 10. Logarithmic plots of § against s for the luminosity sub-samples defined in Table ^. The inset shows a linear plot of § against 
s. The error bars are derived from the scatter between ten realizations of the ACDM mock galaxy catalogues. For clarity, error bars are 
plotted only for the bright sub-sample. 



lar split into low, medium and high 0[II] EW samples gives 
so values of 8.6, 4.9 and 4.1 /t _1 Mpc respectively. The low 
EW samples are early type galaxies, which do not appear in 
IRAS samples. The medium and high EW samples are ac- 
tively star-forming galaxies, which are more similar to IRAS 
samples. The clustering amplitude of the medium and high 
EW galaxies is very similar to that of the IRAS galaxies, 
suggesting that they do indeed trace a similar population. 
Furthermore, the change in so between the medium and high 
SAPM emission-line galaxies is very similar to the change 
between cold and hot PSCz galaxies as seen in Figure |^. 

The luminosity dependence is a more complex effect, 
because the FIR luminosity of a galaxy depends on both the 
mass of dust, and the recent star-formation rate. If the dust 
mass is correlated with the halo mass, then more luminous 
galaxies will, on average, be in more massive halos, which 
are expected to have a higher clustering amplitude. On the 
other hand, a high star-formation rate will make a galaxy 
more luminous, and galaxies with a high star-formation rate 
tend to be less clustered. It is not obvious which will be the 
dominant effect. 

Given the correlation between colour and luminosity 
seen in Figure ^| we can infer the luminosity dependence 
from the colour dependence seen in Figure 0. From Fig- 
ure | we find log 10 (^) = -0.117 log 10 (L 6 o) + 1-427, and 
so together with equation ^ we would expect log 10 (so) = 
—0.028 log 10 (I/6o) + 0.95. This relation is plotted as the 
dashed line on Figure [l2| and it can be seen that it is 
consistent with the data points. The luminosity variation 
seen in optical galaxy samples is approximately of the form 



log 10 (s ) = 4.37(0.7 + 0.3-^) (Benoist et al. 1996). The dot- 
ted line shows this with L* = 3.6 x 1O 9 L (Springel & White 
1998). This is clearly inconsistent with the trend seen in the 
PSCz sample. 

The small decrease in clustering that we find for more 
luminous galaxies is consistent with the colour dependence 
effect alone, with no evidence for an increase in clustering 
for more luminous galaxies. If the standard picture of high- 
peak biasing is correct, we conclude that there is a very weak 
correlation between the FIR luminosity of a galaxy and the 
mass of it's dark-matter halo. However, models which allow 
a variable number of galaxies per dark halo, and include 
scatter about the mean mass-luminosity relation lead to a 
very weak luminosity dependence of clustering (Somerville 
et al. 1999). 

Finally we note that our present analysis has been re- 
stricted to the redshift space correlation functions. This 
means that any differences in the velocity distributions of 
the different sub-samples of galaxies will affect the ampli- 
tude of £(s), particularly on small scales where the velocity 
separations are comparable to the random peculiar veloci- 
ties. We consider these effects in a future paper. 
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Figure 11. Clustering length so versus volume radius for the volume limited sub-samples for a fixed 7 = 1.30. The error bars show the 
overall uncertainties including both Poisson errors in the pair counts, and an estimate of the cosmic variance as described in Section 3.2. 
For volumes with radius R < 100 h~ -"-Mpc we also estimated the uncertainties from the standard deviation between the ten realizations 
of the ACDM mock catalogues and found good agreement with the analytic estimates. 
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Figure 12. Log plot of clustering length so versus mean absolute luminosity for the volume limited sub-samples (filled circles), the bright 
and faint sub-samples (open circles) and the whole sample (open square). The dashed line shows the trend expected from the colour 
dependence seen in Figure ^ and the luminosity-colour correlation in Figure ^ The dotted line shows log 10 (so) = 4.37(0.7 + 0.3-^), 
which approximates the variation seen in optical galaxy samples. 
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